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ABSTRACT 



Type la supernovae (SN la) are the most important standard candles for measuring the expansion history of the universe. The ther- 
^ . monuclear explosion of a white dwarf can explain their observed properties, but neither the progenitor systems nor any stellar remnants 

O^ ' have been conclusively identified. Underluminous SN la have been proposed to originate from a so-called double-detonation of a white 

lO ' dwarf. After a critical amount of helium is deposited on the surface through accretion from a close companion, the helium is ignited 

^^ , causing a detonation wave that triggers the explosion of the white dwarf itself. We have discovered both shallow transits and eclipses 

^^ . in the tight binary system CD-30° 1 1223 composed of a carbon/oxygen white dwarf and a hot helium star, allowing us to determine 

its component masses and fundamental parameters. In the future the system will transfer mass from the helium star to the white dwarf. 

Modelling this process we find that the detonation in the accreted helium layer is sufficiently strong to trigger the explosion of the 
^Tjf , core. The helium star will then be ejected at so large a velocity that it will escape the Galaxy. The predicted properties of this remnant 

are an excellent match to the so-called hypervelocity star US 708, a hot, helium-rich star moving at more than 750 km s"', sufficient 

to leave the Galaxy. The identification of both progenitor and remnant provides a consistent picture of the formation and evolution of 

underluminous type la supernovae. 



Key words, binaries: spectroscopic - subdwarfs - supernovae 

1. Introduction nels, but not conclusively identified yet. Although most SNIa 

form a homogeneous class, about one third of them differ signif- 

The search for the progenitors of SN la is ongoing, but the obser- jcantly in their luminosities and other observational properties 

vational evidence remains inconclusive. In the standard single- ^nd their proper classification is crucial when usi ng suc h events 

degenerate scenarios mass is transfeiTed in a stable way by either ^s standard candles for cosmology (Wang & Han'Sora. 

a main sequence star or a red giant to a white dwarf (WD) com- Underluminous SN la have been proposed to originate from 

panion. In the double-degenerate scenario a close binary consist- ^ so-called double-detonation of a white dwarf. After a critical 

ing of two white dwarfs shrinks due to angular momentum lost ^j^ount of helium is deposited on the surface through accretion 

by the emission of gravitational waves and eventually merges, f^.^^ ^ ^.j^.g companion, the helium is ignited causing a deto- 

Possible progenitor systems have been proposed for both chan- ^^^^^^ ^^^g j^at triggers the explosion of the white dwaif itself 

even if its mass is significantly lower than the Chandrasekhar 

Send offprint requests to: S. Geier, limit (Nomoto[T982l Woosley et al. |1986l l. Hydrodynamic simu- 

e-mail: sgeier@eso . org lations predict the explosion of a CO- WD with a minimum mass 
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of only ~ 0.8 Mq as underluminous SN la triggered by the ig- 
nition of an He-shell of ~ 0.1 Mq (Fink et al. '2010^ in this so- 
called double-detonation scenario. He-stars have already been 
proposed as possible donors for the single-degenerate scenario 
(Yoon & Langer l2003l Wang et al. |2009a| |20^ conveniently 
explaining the lack of hydrogen in the spectra of SN la. Recent 
studies indicate that this scenario might also be consistent with 
the lack of helium in standard SN la spectra as long as the ac- 
creted He-layer is thin (Sim et al. 120 101 Kromer et al. l2010) . 

In the course of the MUCHFUSS project (Geier et al. lSOTlT l. 
which aims at finding hot subdwarf binary systems with mas- 
sive companions, we have discovered a possible progenitor for 
such a supernova consisting of a hot subdwarf B star (sdB) and a 
white dwarf in an extremely compact binary (Heber et al. 120131 
Geier et al. I2012I I. This system, CD-30°11223, was also inde- 
pendently discovered by Vennes et al. ( I2012I I. 

Hot subdwarf stars are evolved, core helium-burning objects. 
About half of the sdB stars reside in close binaries with periods 
rangi ng from ~ 0.1 d to ~ 30d (Maxted et al. 1200 II Napiwotzki 
et al, I2004al l. The sdB is the core of a former red giant star that 
has been stripped off almost all of its hydrogen envelope through 
interaction with a companion star The mass of the emerging sdB 
star is constrained to about half a solar mass in order to allow 
central helium burning. After the helium-burning phase the sdB 
star will turn into a white dwarf. 

Because the components' separation in these systems is 
much less than the size of the subdwarf progenitor in its red-giant 
phase, these systems must have experienced a common-envelope 
(CE) and spiral-in phase (Han et al. 20021 12003] and references 
therein). In this scenario, a star evolving to become a red giant 
swallows a nearby companion. Due to friction the core of the gi- 
ant and the more compact companion spiral towards each other 
in a common envelope. The orbital energy lost during this pro- 
cess is deposited in the envelope until it is eventually ejected 
leaving a close binary system as remnant. 

Although most of the close companions to sdB stars are 
low-mass main sequence stars, brown dwarfs or low-mass WDs 
(~ 0.5 Mq), more massive compact companions like WDs, neu- 
tron stars or black holes have been either observed or predicted 
by theory (Geier et al. 20071 WWi . The short-period sdB -i- WD 
binary KPD 1930-1-2752 is regarded as progenitor candidate for 
an SN la (Maxted et al. lIOOOl Geier et al. l2U07] l. 

Here we report the discovery of both shallow transits and 
eclipses in the tight binary system CD-30° 1 1223 composed of a 
carbon/oxygen white dwarf and a hot helium star, allowing us to 
determine its component masses and fundamental parameters. 
This system turns out to be an excellent progenitor candidate 
for the double-detonation SN la scenario and can be linked to 
the hypervelocity subdwarf US 708, the likely donor remnant of 
such an event. 

2. Observations 

CD-30°11223 (ffiooo = 14*'11'"16?2, ^2000 = -30°53'03", 
niy - 12.3 mag) was selected as UV-excess object and spec- 
troscopically identified to be an sdB star (Vennes et al. 120111 
Nemeth et al. I2012I I. We selected this star as a bright backup 
target for our MUCHFUSS follow-up campaign. Due to un- 
favourable observing conditions, which prevented us from ob- 
serving our main targets, two medium resolution spectra (R ~ 
2200, /I = 4450 - 5 110 A) were taken consecutively with the 
EFOSC2 spectrograph mounted at the ESONTT on June 10, 
2012. The radial velocity shift between those two spectra turned 
out to be as high as 600 km s"^ 




3800 



3900 4000 4100 4200 

Wavelength [Angstroms] 



4300 




60 

T 40 

E 20 

" 

cr -20 
< 

-40 



H 1 1 1 1 1 1 1 1 h 



H 1 1 H 






0.5 



1 
phase 



1.5 



Fig. 1. Upper panel: Spectrum of CD-30°1 1223 coadded from 
175 RV-corrected spectra taken with SOAR/Goodman. The hy- 
drogen Balmer series is clearly visible as well as a prominent 
helium line at 4026 A. The other features are rotationally broad- 
ened metal lines. Model spectra were matched to the hydrogen 
and helium lines to determine the atmospheric parameters of 
the hot subdwarf star. Middle panel: Trailed spectra taken with 
SOAR/Goodman and showing the short-period sinusoidal varia- 
tions of the Doppler-shifted spectral lines caused by the motion 
of the visible sdB star. The close and compact white-dwarf com- 
panion of the subdwarf star is not visible in this dataset. Lower 
panel: Radial velocity curve of CD -30° 1 1223 derived from 105 
spectra taken with WHT/ISIS plotted twice against orbital phase 
for better visualisation. 



First spectroscopic follow-up data was obtained with the 
grating spectrograph mounted on the SAA0-1.9m telescope on 
July 2, 2012. The RV-curve derived from 18 single spectra 
confirmed the short orbital period of 0.0498 d and a high RV- 
semiamplitude (K - 370 + 14kms"'). In order to improve 
the orbital solution and minimize the effect of orbital smear- 
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Fig. 2. Fit of synthetic LTE models to the hydrogen Balmer lines Fig. 3. Fit of synthetic LTE models to the helium lines (see 
of a coadded ISIS spectrum. The normalized fluxes of the single Fig|2]i. 
lines are shifted for better visualisation. 



ing, we took another 105 spectra with the ISIS spectrograph 
(R ~ 4000,^ = 3440 - 5270A,rexp = 2min) mounted at the 
WHT during a dedicated MUCHFUSS follow-up run from July 
9 to 12, 2012. Another 175 spectra were taken with the Goodman 
spectrograph mounted at the SOAR telescope (R ~ 7700, A - 
3700 - 4400 A, Texp = 1 min) on July 16, 2012. 

CD-30° 11223 was observ ed by the SuperWASP planetary 
transit survey (Pollacco et al. '2006'l. The light curve contains 
23 678 measurements taken from May 4, 2006 to August 1, 
2011. 

3.6 hr of time-series photometry in the V-band (Texp = 3 s) 
were taken with SOAR/Goodman on July 6, 2012 under photo- 
metric conditions. The light curve was extracted using an aper- 
ture that minimizes the standard deviation of the two comparison 
stars used divided by each other at low airmass and a flat-field 
correction has been applied. The combination of short expo- 
sure times and bright stars make scintillation the dominant noise 
source for our photometry. Another source of noise is likely to 
be caused by a small scatter in effective integration time of the 
order of a few ms. 



3. Orbital and atmospheric parameters 

The light curve shows variations caused by the ellipsoidal de- 
formation of the sdB primary, which is triggered by the tidal 
influence of the compact companion, as well as Doppler boost- 
ing, caused by the extreme orbital motion of the sdB (Shakura 
& Postnov 119871 see also Bloemen et al. 120111 and references 
therein). The ephemeris has been derived from the SWASP data 
based on fitting a harmonic series. Due to the timebase of more 
than five years the derived orbital period of 0.0489790724 + 
0.0000000018 d is very accurate and perfectly consistent with 
the independent determination (P = 0.04897906+0.00000004 d) 
by Vennes et al. (120 12l . 

The radial velocities were measured by fitting a set of math- 
ematical functions to the hydrogen Balmer lines as well as he- 
lium lines using the FITSB2 routine (Napiwotzki et al. l2004b] l. 
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Fig. 4. Teff - \ogg diagram. Evolutionary tracks (solar metal- 
licity) of core helium-burning star with masses of 0.45 Mq 
(dotted Unes), 0.50 Mq (short-dashed fines) and 0.55 Mq (long- 
dashed lines) are plotted for different hydrogen envelope masses 
(0.000 Mo, 0.001 Mo, 0.005 Mq from the lower left to the upper 
right).* The diamond marks CD-30°1 1223. 



Three functions (Gaussians, Lorentzians and polynomials) are 
used to match the continuum, the line and the line core, respec- 
tively and mimic the typical Voigt profile of spectral lines. The 
profiles are fitted to all suitable lines simultaneously using x^- 
minimization and the RV shift with respect to the rest wave- 
lengths is measured. Assuming circular orbits sine curves were 
fitted to the RV data points in fine steps over a range of test pe- 
riods. The two datasets obtained with ISIS and Goodman are 
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Fig. 5. Upper panel: V-band light curve of CD -30° 11223 taken with SOAR/Goodman (green) with superimposed model (red) 
plotted twice against orbital phase for better visualisation. The dashed red curve marks the same model without transits and eclipses. 
The sinusoidal variation is caused by the ellipsoidal deformation of the hot subdwarf due to the tidal influence of the compact white 
dwarf. The diff'erence in the maxima between phase 0.25 and 0.75 originates from the relativistic Doppler boosting eff'ect, which is 
usually not detectable with ground-based telescopes. Lower panels: Close-up on the transit of the WD in front of the sdB (left). It is 
even possible to detect the eclipse of the WD by the sdB (right). 



treated separately to investigate systematic errors. Details about 
the analysis method and error estimation are given in (Geier et 
al.|20H). The derived orbital parameters from the ISIS dataset 
{K = 378.6 + l.Okms-i, y = 17.6 + 0.7 km s"', see Fig. [T] 
lower panel) and the Goodman dataset {K - 374.5 ±1.1 km s ', 
y — 21.3+0.8kms ', see Fig.lT] upper panel) are consistent tak- 
ing into account that systematic uncertainties are usually some- 
what higher than the statistical 1 cr errors given here. The devi- 
ation in system velocity is most likely caused by a slight sys- 
tematic zero-point shift between the two instruments. For fur- 
ther analysis we used the average values. Those values are in 
reasonable agreement with the results {K — 386.9 + 1.9kms ', 
y = 31.5 + 1.3 km s') of Vennes et al. (120121 1. although some- 
what discrepant. 

The atmospheric parameters efifective temperature T^s, sur- 
face gravity log g, helium abundance and projected rotational 
velocity were determined by fitting simultaneously the observed 
hydrogen and helium lines of the single spectra with metal- 
line-blanketed LTE model spectra (Heber et al. 120001) as de- 
scribed in Geier et al. (I2007I I. No significant variations of the 
parameters with orbital phase have been detected. Average val- 
ues and standard deviations have been calculated for the ISIS 
(T^eff = 28800 + 200 K, log ^ = 5.67 +0.03, logy = -1.50 + 0.07, 
Vrotsin/ = 180 + 8kms"') and Goodman datasets (reff - 
29600 + 300 K, \ogg = 5.65 ± 0.05, logy = -1.46 + 0.14, 
Vjotsin/ = 174 + 12kms"'), separately. We adopt the average 
values from both datasets for further analysis. The final helium 
abundance is taken from the ISIS data, because of the higher 
number of He-lines in the spectral range. 



The derived parameters are consistent with literature values 
within the uncertainties (Vennes et al. '2011', Nemeth et al.l 2012l l. 
More detailed information about the systematic errors of this 
method can be found in Geier et al. (2007 201 1). Table 1 shows 
the orbital and atmospheric parameters. Fig. |4]the position of 
CD-30° 1 1223 in the T^s - log g diagram. 



4. Light curve analysis 

The light curve obtained with SOAR/Goodman was analysed by 
fitting models calculated with the lcurve code written by TRM 
(see Fig. |5] Copperwheat et al. 1201 Oi l. The code uses grids of 
points modelling the two stars and takes into account limb dark- 
ening, gravity darkening, mutual illumination effects, Doppler 
boosting and gravitational lensing. Since the masses and radii 
of both components are strongly correlated, those parameters 
have been constrained using Markov chain Monte Carlo simu- 
lations. A detailed description of the analysis method is given by 
Bloemen et al. (120111) . 

In order to determine masses and radii of both the sdB and 
the WD, we used two different prior constraints. In each case 
the K derived from spectroscopy is used. First, we assumed tidal 
synchronisation of the sdB primary (solution 1), which is a rea- 
sonable assumption given the short orbital period of the system. 
The Vrot sin / derived from spectroscopy is measured making the 
simplified assumption of a spherical, linear limb darkened star 
(limb darkening coefficient 0.3). In order to take into account the 
additional eff'ects of limb darkening and gravitational darkening 
of the Roche-distorted star we calculated a correction factor of 
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0.963 by comparing the slightly different Une profiles calculated 
under both assumptions (Claret & Bloemen'201D. This correc- 
tion was applied to the Viot sin / before deriving the binary param- 
eters. The best fit is achieved for an sdB mass of 0.47 + 0.03 Mq 
and a WD mass of 0.74 + 0.02 Mq. However, we found that the 
radius of the WD is about 10% smaller than predicted by the 
zero-temperature mass-radius relation for WDs, which provides 
a lower limit for the WD radius (see Fig. |6] second panel from 
above, Verbunt & Rappaport |1988l l. 

In order to explore the influence of this discrepancy, we im- 
posed the restriction that the white dwarf is within 2% of the 
M-7?-relation and allowed for deviations from corotation (so- 
lution 2). We determined an estimate for the temperature of 
the WD from our light curve analysis. Since we see a signif- 
icant feature when the sdB occults the WD, we can derive a 
black-body temperature of 24700 + 1200 K for the WD, which 
leads to a radius about 5% higher than expected from the zero- 
temperature relation. We adopt this more realistic value for our 
analysis. In this case the derived masses are somewhat higher 
(MsdB = 0.54 + 0.02 Mo, MwD = 0.79 + 0.01 Mq). 

Although both solutions are consistent within their uncer- 
tainties we refrain from favouring one over the other In order 
to calculate the kinematics and the further evolution of the sys- 
tem, we adopt the average values for the component masses 
(MsdB = 0.51 Mo, MwD = 0.76 Mo, see Fig. IS. 

Comparing the derived sdB masses with evolutionary tracks 
for core helium-burning objects (Fig. |4}, it can be seen that 
the appropriate tracks are consistent with the position of 
CD-30° 11223 in the Teff - \ogg diagram. Furthermore, the ef- 
fective temperature and surface gravity of the star tell us that the 
sdB has just recently been formed and started the core helium- 
burning phase, which typically lasts for about 100 Myr. 

5. Gravitational wave radiation 

Because of its short orbital period, CD -30° 1 1223 is expected to 
be a strong source of gravitational waves. We therefore calcu- 
lated the current gravitational wave emission of CD-30°1 1223. 
The gravitational wave strain amplitude h scales with the masses 
of both binary components, the binary inclination, the orbital pe- 
riod and the distance of the system. 

We calculate it as described in Roelofs et al. (I2007I I to be as 
high as log/i = -21.5 + 0.3. CD-30°11223 should therefore 
be one of the strongest gravitational wave sources detectable 
with missions like NGO/eLISA (Kilic et al. 120 12t Nelemans 
12009 ). It sticks out even further, because due to the presence 
of eclipses its binary parameters are determined to very high ac- 
curacy. Therefore this system can be used as verification source 
for upcoming space missions. 

No period change due to the orbital shrinkage caused by the 
emission of gravitational wave radiation has been detected in the 
SWASPdata(P= 1.01 x lO^'^ + 3.38 x lO^'^ ss^^). This non- 
detection is consistent with the theoretically expected value of 
P ~ 6 X 10 '"* s s '. However, within only a few more years the 
orbital shrinkage should become detectable. 



6. Binary evolution calculations 

More interesting than the present state of this system is its future 
evolution, which can now be studied in detail using theoretical 
models. Employing Eggleton's stellar evolution code (Eggleton 
1197 1 1972 1973 ), we calculate the evolution of the sdB star and 
its WD companion. The code has been updated with the latest in- 
put physics over the past four decades (Han et al. il994i Pols et 



al. [19951 [T998] l. Roche lobe ove rflow (RLOF) is treated within 
the code described by Han et al. (2000). We set the ratio of mix- 
ing length to local pressure scale height, a - l/Hp, to be 2.0. We 
assume that the binary model starts with a 0.51 Mo He star and a 
0.76 Mo CO WD having a 0.049 d orbit period, similar to the ini- 
tial model of the sdB star and its WD companion. Additionally, 
orbital angular momentum loss due to gravitational wave radi- 
ation is included by adopting a standard formula presented by 
Landau & Lifshitz (fTgTTT l. 



d In Jgr 32G^ MwdM2(Mwd + M2) 



dt 



5c5 



(1) 



where G, c, Mwd and Mi are the gravitational constant, vacuum 
speed of light, the mass of the accreting WD and the mass of the 
companion sdB star, respectively. 

In the He double-detonation model, if the mass-accretion 
rate is higher than 4 x lO^^Moyr ', the WD can increase its 
mass (Woosley et al. '1986'; Wang et al. ■2009a'). However, for 
low mass-accretion rates (1 x lO^^Moyr"' < IM2I < 4 x 
10"** Moyr"'), compressional heating at the base of the accreted 
He layer plays no significant role, and a layer of unbumed He 
can be accumulated on the surface of the WD. For the case with 
a slower mass-accretion rate (IM2I < 1 x lO^^Moyr"'), the He 
flash is strong enough to form a He detonation but too weak to 
initiate a carbon detonation, resulting in that only a single He 
detonation wave propagates outward (Nomoto 1 9821) . 

We assume that if such a CO WD accumulating He enters 
this 'low' accretion rate regime and accumulates 0.1 Mo of He 
on its surface, a detonation is initiated at the base of the He shell 
layer (Fink et al. 120101 ). Consequently, a detonation in the core 
of the CO WD is presumed to follow, in which a sub-Mch SN 
la takes place. In this model, only accreting WDs with a total 
mass (CO core 0.8 Mo + heHum shefl 0.1 Mq) >0.9 Mo are 
considered to result in potential sub-Mch SNe la, since the CO- 
core with lower mass may not detonate and it is unlikely to pro- 
duce enough radioactive nickel observed in SNe la (Kromer et 
al. 2010). 

The binary system starts with (M^, M^^, log(PVd)) = 
(0.510, 0.760, -1.310), where M[, M^^ are the initial masses 
of the He star star and of the CO WD in solar mass, and P' is 
the initial orbital period in days. Fig. |7] (left panel) shows the 
evolutionary track of the He star and the evolution of the orbital 
period. Fig. Q (right panel) displays the mass-transfer rate and 
the mass of the WD envelope varying with time after the He star 
fills its Roche lobe. 

Due to the short initial orbital period (0.049 d) of the sys- 
tem, angular momentum loss induced by gravitational wave ra- 
diation is large. This leads to the rapid shrinking of the or- 
bital separation. After about 36 million years, the He star be- 
gins to fill its Roche lobe while it is still in the core helium- 
burning stage. The mass-transfer rate is stable and at a low rate 
between 1.6 x 10"**Moyr"' and 2.2 x 10"**Moyr"', resulting 
in the formation of a He shell on the surface of the CO WD. 
After about 6 million years, the mass of the He shell increases to 
~ 0.1 Mq in which a double-detonation may happen at the base 
of the He shell layer. At this moment, the mass of the He star is 
Mf* = 0.41 Mq and the orbital period is log(P™/d) = -1.72 
(P = 0.019 d). 



7. Hypervelocity sdO as donor remnant 

Theoretical predictions about whether or not a progenitor candi- 
date will explode as SN la are useful, but in general difficult to 
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Fig. 6. Parameters of the CD-30° 1 1223 system plotted against the mass of the WD companion with contours marking the 1 cr (red), 
2 o" (light blue), and 3 cr (dark blue) levels of confidence. In the left column tidal synchronisation of the sdB primary is assumend, 
while in the middle column the effect of microlensing has been taken into account in addition. In the right column the WD mass- 
radius relation has been used as additional constraint instead. The solid lines in the mass-radius plots (second panels from above) 
mark the zero-temperature mass-radius relation for WDs, the dashed lines more realistic models with 5% inflation. Note that both 
the v'lot sin / and the log g derived from the light curve are in agreement with the values derived from spectroscopy. 



test. Usually the theoretically predicted SN rates are compared to 
the observed ones, but these comparisons are often hampered by 
selection effects. A more direct proof would be the identification 
of the remnant objects. We therefore follow the future evolution 
of CD-30°11223. At the end of the He-accretion phase and just 
before the SN event, the orbital period of the binary is predicted 
to have shrunk to 0.019 d due to the further loss of orbital energy 
through the emission of gravitational waves. The sdB primary 
lost a fair amount of mass (~ 0.1 Mq), which was transferred 
to the WD companion. The orbital velocity of the sdB will be 
about 600 km s"' and therefore close to the Galactic escape ve- 
locity. As soon as the WD is disrupted, the sdB will be ejected. 
Depending on the ejection direction of such an object relative to 
its trajectory around the Galactic centre, the Galactic rest frame 
velocity could be even higher by up to 240 km s"'. In this case 
the remnant star will leave the Galaxy. 

Such so-called hypervelocity stars have indeed been dis- 
covered (Brown et al. 120051 Hirsch et al. 2005 Edelmann 
et al. 120051) . However, all but one of the known 22 objects 
are intermediate-mass main-sequence star. This enigmatic star 



(US 708) has been classified as helium-rich hot subdwarf trav- 
elling at a Galactic rest frame velocity of at least 750km s"' 
(Hirsch et al. 2005), which matches the predicted ejection ve- 
locity of CD-30°11223 very well. It was proposed that this 
star might be the ejected He-donor after the WD companion ex- 
ploded as SN la (Jusfliam et al. |2009l Wang & Han |2U09] l. 

In this scenario, the compact binary CD-30° 11223 and the 
hypervelocity star US 708 represent two different stages of an 
evolutionary sequence linked by a SNIa explosion. The exis- 
tence of objects like US 708 thus provides evidence that binaries 
like CD-30°1 1223 are viable SN la progenitor candidates. 



8. Age of the binary system 

The analysis of our data also allows us to constrain the initial 
component masses and the age of the binary. Furthermore we 
can constrain both its past and future trajectory. 
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Table 1. Parameters of the CD-30° 11223 system 



Visual magnitude'' 


my 


[mag] 


12.342 + 0.003 


Proper motion* 


yu„ cos 6 


[mas/yr] 


9.5 ±2.2 




tl6 


[mas/yr] 


-5.6 ±2.2 


Atmospheric parameters of the subdwarf 


Effective temperature 


Teff 


[K] 


29 200 ±400 


Surface gravity 


log,? 




5.66 ±0.05 


Helium abundance 


l0g.V 




-1.50 ±0.07 


Projected rotational velocity 


v,ot sin ;■ 


[kms-'] 


177 ± 10 


Orbital parameters 




To 


[BJD UTC] 


2455113.205908 ±0.000363 


Orbital period 


P 


[d] 


0.0489790724 ± 0.0000000018 


RV semi-amplitude 


K 


[Ions"'] 


376.6 ± 1.0 


System velocity 


y 


[Ions-'] 


19.5 ± 2.0 


Binary mass function 


f(M) 


[Mo] 


0.271 ±0.002 


Derived parameters 


Solution 1 


sdB mass 


MsiB 


[Mo] 


0.47 ± 0.03 


sdB radius 


RsdB 


[^o] 


0.169 ±0.005 


WD mass 


MwD 


[Mo] 


0.74 ± 0.02 


WD radius 


^WD 


[Re] 


0.0100 ±0.0004 


Orbital inclination 


i 


n 


83.8 ± 0.6 


Separation 


a 


[Re] 


0.599 ± 0.009 


Mass ratio 


9 




0.63 ± 0.02 


Solution 2 


sdB mass 


MsdB 


[Me] 


0.54 ± 0.02 


sdB radius 


^sdB 


[Re] 


0.179 ±0.003 


WD mass 


MwD 


[Me] 


0.79 ± 0.01 


WD radius 


^WD 


[Re] 


0.0106 ± 0.0002 


Orbital inclination 


i 


[°] 


82.9 ± 0.4 


Separation 


a 


[^o] 


0.619 ± 0.005 


Mass ratio 


1 




0.68 ± 0.01 










ide is taken from Vennes et al. (1201 2l 




en from Roeser et al. pOlO). 









8.1. Kinematic analysis 

Using a standard Galactic gravitational potential with a Sun- 
Galactic centre distance of 8.4 kpc and a local standard of 
rest circular motion of 242 kms"' (see Model I in Irrgang 
et al. I2013I I. we computed the past and future trajectory of 
CD-30° 11223 (see Fig. |8]l. The orbit shows the typical char- 
acteristics of the local thin disc population, i.e., almost circular 
motion around the Galactic centre and small oscillations in di- 
rection perpendicular to the Galactic disc. The heliocentric dis- 
tance to the star increases during the next 42 Myr until the super- 
nova is predicted to explode from its current value of 364 + 3 1 pc 
to about 1920+ 160 pc. 

CD-30° 11223 is by far the closest known SNIa progenitor 
with respect to Earth. The explosion will take place in a direc- 
tion of the sky close to the current positions of the constellations 
Ara and Norma. Adopting an absolute visual magnitude of up to 
- 19 mag for the SN la, the apparent magnitude seen from Eaith 
might be as high as ~ -7.6 mag or about as blight as SN 1006, 
the brightest stellar event in recorded history so far (Winkler et 
al. 12003,) . 



8.2. Binary formation scenario 

CD-30° 11223 is the closest sdB binary known so far and the 
mass of its WD companion is higher than the average mass of 
CO WDs (~ 0.6 Mq). In order to explore the formation of this 
exceptional system, we performed a binary population synthesis 
study in a similar way as described by Han et al. (2002, 2003). 
For given WD masses ranging from 0.6 Mq to 1 .4 Mq, an initial 
set of 10^ WD-i-MS binaries was generated. For the main se- 
quence stars the initial mass function of Salpeter was used. The 
orbital peiiod distiibution was assumed to be flat in log a. The 
binaries have been evolved through the common envelope phase 
for different values of the CE-efficiency parameters acE, which 
is the fraction of the available orbital energy used to eject the 
envelope, and Cth, the contributed fraction of internal energy. 

In the standard scenario, which is called the 2nd CE channel, 
the progenitor of the sdB is a main sequence star of about solar 
mass and the common envelope is ejected right at the tip of the 
first giant branch (FGB, Han et al. 2002 ). However, this channel 
is not feasible to form binaries as close as CD-30°11223, as 
shown in Fig. |9](upper panel). The envelope at the tip of the FGB 
has a very low binding energy and can be ejected easily in the 
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Fig. 7. Results of binary evolution calculations with initial masses of the two components and the orbital period similar to the 
sdB+WD binary system. Left panel: The evolutionary track of the He star is shown as solid curve and the evolution of the orbital 
period as dash-dotted curve. Right panel: The solid and dash-dotted curves show the mass-transfer rate and the mass of the WD 
envelope (He shell) varying with time after the He star fills its Roche lobe, respectively. Dotted vertical lines in both panels indicate 
the position where the double-detonation may happen (the mass of the He shell increases to ~ 0. 1 Mq). The initial binary parameters 
and the parameters at the moment of the SN explosion are also given in the two panels. 



0.25 kpc 




Fig. 8. Three-dimensional trajectory of CD-30° 11223 in a 
Cartesian Galactic coordinate system with the z-axis pointing 
to the North Galactic pole. Current positions of CD-30° 11223 
(blue •), Sun (black 0), and Galactic centre (black +) are 
marked. The approximate point in time of the supernova explo- 
sion is symbolized by the red asterisk, while the arrow marks the 
position of the Sun at that time. Solid lines indicate the future 
42Myr, dashed lines the past 150 Myr CD-30°11223's kine- 
matic properties are obviously those of the local thin disc popu- 
lation. 



following CE. Thus, the orbital shrinkage during CE evolution 
is not significant and the produced sdB+WD system generally 
has an orbital period much longer than that of CD-30° 11223. 



Only for a very small value of acE — 0.3, which is very unlikely, 
some binaries reach the margin of 0.05 d. Indeed, the median 
period of the observed sdB binaries is as high as ~ 0.6 d (Geier 
etal. l20m) . 

However, an sdB+WD binary can also be formed when the 
main-sequence progenitor of the subdwarf has an initial mass 
larger than 2 M© and fills its Roche lobe during the Hertzsprung 
Gap or at the base of the FGB. In this case, the envelope is more 
tightly bound and the orbital shrinkage required to eject the CE 
becomes higher In Fig. |9](lower panel) the orbital period distri- 
bution is shown for this scenario when acE = Cth - 0.75, similar 
to the best fitting model of Han et al. (i2003l) . As seen in the fig- 
ure, short orbital periods just as in the case of CD -30° 1 1223 are 
expected. 

Additional to the orbital period distribution, we also inves- 
tigated the distribution of sdB masses formed via this channel. 
While the standard CE-scenario predicts a mass distribution with 
a sharp peak at 0.47 M©, the sdB masses from more massive 
main-sequence stars (i.e. > 2 M^^) show a significant scatter for 
higher values of acE and even more so, if we allow for a con- 
tribution of thermal energy in the CE-process by increasing the 
parameter ath- The sdB mass for this channel largely depends on 
the mass of the progenitor and can range from 0.3 to 1 .0 M© (see 
Fig. [Toll. This is consistent with the sdB mass of up to 0.54 Mq 
determined in the case of CD-30° 11223. 

We therefore conclude that CD-30° 11223 was most hkely 
formed via CE-ejection of a main sequence star with a mass 
larger than 2 Mq, which means that it originated from a young 
stellar population. 

8.3. White dwarf cooling age and progenitor masses 

We derived a mass of 0.76 Mq for the WD companion based on 
observations. Using an initial-to-final mass relation for isolated 
WDs the mass of the progenitor should have ranged from 3 to 
4Mq (see Fig. 7 in Kovetz et al. (2009) and references therein). 
Binary evolution is expected to lower the mass of the final WD 
and this progenitor mass estimate therefore has to be considered 
as lower limit. 
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Fig. 9. Top panel: Orbital period distribution of sdB+WD bi- 
naries from the 2nd CE channel. Different colours mark dif- 
ferent WD masses (black 0.6 Mq, red 0.8 Mq, green 1.0 Mq, 
blue 1.2 Mq, purple 1.4 M©). The vertical dashed line marks 
P = 0.05 d. The contribution of thermal energy has not been in- 
cluded (ffth - 0). The solid lines are for uqe = 0.3 and the dotted 
ones for acE = 1 0. Bottom panel: Orbital period distribution of 
sdB+WD binaries from the CE channel involving more massive 
main sequence stars (i.e. > 2 Mq). The values of acE and a^i are 
indicated in the figure (see Fig.|9]l. 



Assuming a lifetime on the main sequence tms - lO'^yr x 
[M/Mq]^^-^ the progenitor lived for a maximum of 640 Myr. We 
constrain the temperature of the WD in our light curve analysis 
to be ~ 25000 K. Therefore its cooHng age is ~ 30 - 40My£]. 
The lifetime of the sdB on the extreme horizontal branch is of the 
same order as the one of the WD and we therefore derive a sim- 
ilar progenitor mass of more than 3 Mq consistent with the disc 
kinematics and the constraints from binary formation scenarios. 



' Taken from cooling models for mixed C/0 com- 
position and for standard thick H and He layers, 
|http://www.astro.umontreal.ca/~bergeron/CoolingModels/| Fontaine et 
al. 120011 



log P(days) 

Fig. 10. sdB mass is plotted over the orbital period assuming the 
mass of the WD to be 1.2 Mq. The main-sequence progenitors 
of the sdBs have initial masses larger than 2 Mq. The values of 
acE and ath are indicated in the figure. 



9. Conclusions 

Systems like CD-30° 1 1223 are young, which is consistent with 
the non-detection of objects with such high RV-shifts in the 
course of the MUCHFUSS project so far (Geier et al. 120111 
12011) . Most targets of this survey are faint subdwarfs located 
in the old halo population. The sdB-donor double-detonation 
channel is therefore predicted to occur in young stellar popu- 
lations and contribute to the SNIa population with short delay 
time (Ruiter et al. 2009 ). 

Given that systems like CD-30° 11223 are progenitors of 
some thermonuclear SN, a rough estimate can be made about 
the rate of such events. CD-30°11223is one out of ~ 100 solved 
sdB binaries (Geier et al. 1201 11 1. About 50% of the known sdB 
stars are in close binary systems. So we can estimate the num- 
ber fraction of systems like CD-30°1 1223 to be about 0.5% of 
the whole population of sdB stars. According to binary evolu- 
tion calculations, the birthrate of such stars in our Galaxy is 
~ 5 X 10"^ yr"' (Han et al. I2003I I. We therefore estimate the 
number of progenitor systems and the resulting SN la rate to be 
~ 2.5 X 10"^ yr"' . This is consistent with the theoretical birthrate 
predicted for the WD-i-He star channel (~ 3 x 10 "^yr', Wang 



et al. 1201 0) But more importantly, it is smaller than the SN la 
birthrate of ~ 3 x 10"^ yr"' and therefore consistent with obser- 
vations (Capellaro & Turato |1997) . 

Although sub-Chandrasekhar scenarios in general have no 
well defined explosion mass, the parameter space for the sdB bi- 
nary progenitors turns out to be quite narrow. According to hy- 
drodynamic simulations the minimum mass of the WD should 
be ~ 0.8 Mq, because carbon burning is not triggered for objects 
of much lower mass (Sim et al. 2012). On the other hand, the 
WD must consist of carbon and oxygen to be able to explode as 
SN la. This limits the mass to values lower than ~ 1.1 Mq, be- 
cause even more massive WDs consist of oxygen, neon and mag- 
nesium and would rather collapse than explode. This mass range 
is further narrowed down by binary evolution calculations. Very 
close sdB+WD systems with companion masses around 0.8 Mq 
are predicted to be formed in much higher numbers than binaries 
with more massive companions. Another important constraint is 
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that the timespan from the binary formation after the CE to the 
SN la explosion must be shorter than the core helium-burning 
lifetime (~ 100 Myr). Otherwise the sdB will turn into a WD be- 
fore helium can be transferred. This restricts the orbital periods 
of possible sdB+WD progenitors to less than ~ 0.07 d. 

The double-detonation scenario with hot subdwarf donor is 
the only proposed SN la scenario where both progenitors and 
remnants have been identified. Analysing a larger sample of 
those objects will allow us to put tight constraints on their prop- 
erties and evolution. 
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